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ABSTRACT 

We construct a starburst model for the hidden starbursts in Arp 220 based on the new Starburst99 
models of Leitherer et al. Comparing these stellar population synthesis models with observations, we 
show that the hidden power source must be due to star formation (as opposed to an AGN) at the 50% 
level or more in order to avoid an ionizing photon excess problem, and this starburst must be young 
(< 7 x 10 7 yr). We derive a current star formation rate of 270M© yr _1 , and an extinction Ay > 30 mag 
for our line of sight to this hidden starburst. 

Subject headings: galaxies: individual (Arp 220) - galaxies: starburst - stars: formation 



1. INTRODUCTION 

Ultraluminous infrared galaxies (ULIGs) are the most 
luminous galaxies in the Universe and their bolometric 
luminosities are comparable to those of quasars (£t>oi > 
lO 12 ^©). The physical origin of the large luminosity of 
ULIGs is unknown since the very centers of these ob- 
jects, which is where the bolometric luminosity comes 
from (Sanders & Mirabel 1996) is very heavily obscured 
by dust and direct observation is not possible. Either an 
enshrouded starburst or an active galactic nucleus (AGN) , 
or both, could generate the observed luminoisites. 

Arp 220 is the nearest and best-studied ULIG 1 Recently 
Genzel et al. (1998) performed a mid-infrared spectroscopy 
program for a sample of ULIGs and demonstrated that 
Arp 220 has an appreciable PAH feature and so is pow- 
ered by a starburst, assuming that the very core where 
-Lboi originates from is not optically thick at 15 /mi. Based 
on detailed radiative transfer modelling, Rowan-Robinson 
(2000 and references therein) argues that galaxies like Arp 
220 must be star-formation powered since they emit most 
of their luminosity longward of 50/Ltm. Additional indi- 
rect evidence that Arp 220 is powered by star formation 
includes 

(i) the lack of an X-ray detection (Iwasawa 1999), even at 
hard wavelengths where Compton scattering is less impor- 
tant; 

(ii) VLBI radio observations, which show many radio su- 
pernovae but no strong point radio source which could 
correspond to the AGN (Smith et al. 1998); 

(iii) a LINER 2 -like optical emission-line spectrum which 
can result from a starburst-driven superwind (Veilleux et 
al. 1995; Taniguchi et al. 1999; Lutz, Veilleux, & Genzel 
1999). This is not, however, unambiguous evidence that 
there is no AGN present. In any case, these measurements 



do not probe the very central regions, which is where most 
of the power is emitted. Rather they probe regions with 
optical depths r < 1 at optical wavelengths. 
On the other hand OH megamaser emission certainly in- 
dicates that the power source is compact (Lonsdale et 
al. 1994, 1998). This would be consistent with an AGN 
power source, but again this measurement is not a unique 
discriminant. The nature of the power source in cold 
ULIGs like Arp 220 has taken on a particular significance 
in recent times with the detection of a population of these 
objects at high redshifts which contribute a very large frac- 
tion of the extragalactic background light at far-infrared 
and submillimctcr wavelengths (Barger, Cowie & Sanders 
1999). 

We now attempt to provide additional constraints on the 
nature of the power source in Arp 220 by investigating the 
so-called ionizing photon deficiency problem. Prestwich 
et al. (1994) performed near-infrared (NIR) spectroscopy 
for a sample of starburst galaxies including Arp 220 and 
found that the ratio of the Br7 luminosity to the infrared 
luminosity of Arp 220 is lower by a factor of ten than 
what is observed in lower-luminosity (but still very lumi- 
nous) starburst galaxies, suggesting a deficit of ionizing 
photons. This could in part be due to absorption of the 
ionizing photons by dust, to an AGN contributing much 
of Lboi but few ionizing photons, or to the detailed shape 
of the stellar initial mass function (IMF). These physical 
process all operate on the global properties of Arp 220 in 
different ways and in this paper we use data from the liter- 
ature together with starburst models and an SED fitting 
method to give constraints on the nature of the hidden 
power source. 



2. MODELS 



1 We use a distance of 74 Mpc, assuming Hq = 75 km s 1 Mpc 1 and Qq = 1 with the recession velocity to the Galactic Standard of Rest, 
Vgsr = 5, 531km s _1 (de Vaucouleurs et al. 1991). 

2 LINER = Low Ionization Nuclear Emission-line Region (Heckman 1980). 
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The nuclear super-star clusters seen in optical (Shaya et 
al. 1994) and near-infrared (Scoville et al. 1998) HST im- 
ages are not the source of the high bolometric luminosity 
of Arp 220 (Shioya, Taniguchi & Trentham 2000), which is 
completely hidden by dust at those wavelengths. It is this 
hidden power source that we attempt to model. We use 
the population synthesis model of Lcitherer et al. (1999, 
Starburst99) to compute bolometric luminosities and ion- 
izing photon fluxes which we compare to the observed val- 
ues derived from far-infrared and radio data respectively 
(where extinction is less of a problem). 

The following assumptions are made: 

1) we assume a model with continuous star formation. 
Note that there is abundant cold gas available as fuel 
(Scoville et al. 1991, 1998; Downes & Solomon 1998); 

2) we consider the three stellar IMFs given in Table 1, 
where we list the IMF slope a, lower-mass cutoff Mi, and 
upper-mass cutoff M u . The values of M\ are higher than 
what is observed in the solar neighborhood, but seems to 
be more appropriate for violent starburst galaxies (e.g., 
Goldader et al. 1997); 

3) we assume solar metallicity, appropriate for the central 
regions of galaxies, which is where most stars generated in 
this kind of dense starburst will end up (e.g. Faber 1973); 

4) we assume that a hypothesized AGN coexists with the 
starburst and contributes a fraction /agn of the bolo- 
metric luminosity. Radio emission from this AGN is as- 
sumed to contribute negligibly to the ionizing photon flux 
inferred from radio recombination lines. This is reason- 
able since we expect an electron column density larger 
than ~ 10 23 cm~ 2 given the existence of very dense gas 
traced by HCN (Solomon, Downes, & Radford 1992) and 
CS (Solomon, Radford, & Downes 1990). Radio emission 
from an AGN is then suppressed because of free-free ab- 
sorption. 

3. RESULTS 

3.1. Initial Mass Function, Age, and AGN fraction 

The total bolometric luminosity of Arp 220 is Lboi — 
1.5 x 10 12 i Q (Sanders et al. 1988). When considering 
the hidden power source, we need to subtract the part 
of the bolometric luminosity which is observed at optical 
and near-infrared wavelengths in the center of the galaxy. 
The K- magnitude of the central 2.5" diameter region of 
Arp 220 is 11.80 mag (Carico et al. 1990). Even with 
Ay = 6.4 mag of extinction (Shier et al. 1996), this is 
only 1 x 10 11 L Q (assuming a bolometric correction of 0.5; 
Leitherer et al. 1999). This value is less than 10 % of the 
total bolometric luminosity of Arp 220 and we neglect it. 

The ionizing photon production rate of Arp 220 is eval- 
uated from the observation of radio recombination lines 
like H92a: Q = 1.3 x 10 55 photons s -1 (Anantharamaiah 
et al. 2000). This is twice as large as the rate derived by 
Zhao et al. (1996) but Anantharamaiah ct al. (2000) have 
wider velocity coverage; note that it is also larger than 
the upper limit derived by Scoville et al. (1991) from the 
free-free radio continuum flux density at 110 GHz, Q < 
7.5 x 10 54 photons s _1 . Again, we have to estimate the 
ionizing photon rate from the regions we do see in if-band 
and subtract it from the above values to derive the pro- 
duction rate for the hidden power source. The reddening- 
corrected (Ay ~ 10 mag; Larkin et al. 1995) Eh-7 flux 



(5.9 x 10~ 18 W m~ 2 ; Goldader et al. 1995) gives an ioniz- 
ing photon rate, 7 x 10 53 photons s _1 for these visible star- 
burst regions. Therefore the production rate for the hid- 
den power source is: Q(hidden) ~ 1.23 x 10 55 photons s _1 . 

Following Scoville et al. (1997), we adopt the 
stellar mass of the nucleus of Arp 220 as 2.5 x 
10 9 M Q . We therefore set the constraints on the stel- 
lar populations of the hidden starbursts as follows: 

8.1 x 10~ 44 L Q (photons s- 1 )- 1 < L bo i/Q(hidden) < 

1.2 x 10~ 43 L (photons s" 1 )" 1 and 400 L Q /M Q < 
Lbo\/M* (hidden) < 600 L Q /M Q . An AGN fraction 
/agn > results in both of these parameters being low- 
ered by a factor (1 — /agn)- 

Now we compare the above numbers to the parameters 
derived for the three models described in Section 2. Fig- 
ure 1 shows the evolutionary locus on a Lboi/Q - -^boi/M* 
plane for each model. The general form of the curves is un- 
derstood as follows. Low-mass stars have smaller L^ \jM^ 
but larger Lboi/Q than high-mass stars. Since the frac- 
tional mass of low-mass stars increases with increasing age, 
as time increases, Lboi/M,, becomes smaller while Lboi/<3 
becomes larger for all the models. 

Let us consider the results for the three IMFs in turn. 

3.1.1. IMF1 

For IMF1, the model overpredicts the bolometric lumi- 
nosity given the observed ionizing photon flux (box A). 
This is the classical ionizing photon deficiency problem 
(e.g., Goldader et al. 1995). This could be solved in a 
number of ways: 

(i) reducing the number of massive stars which generate 
more bolometric luminosity per unit ionizing photon that 
lower mass stars. This could be achieved by either steep- 
ening the slope of the IMF (as in IMF2) or lowering the 
upper mass cutoff (as in IMF3). These adjustments have 
different effects on the properties studied here and are dis- 
cussed in the next two sections. 

(ii) if some ionizing photons are absorbed by dust shortly 
after their production, the ionizing photon rate evaluated 
from radio recombination lines will be too low. The fig- 
ure shows that IMF1 can then be consistent with obser- 
vation if the hidden power source is a pure starburst with 
2 x 10 7 < age/yr < 3 x 10 7 and between 0.3 and 0.6 of the 
ionizing photons that are generated are absorbed by dust 

(iii) if an AGN contributes some of the power, the ionizing 
photon deficiency problem can be reduced. If /agn = 0.5, 
then IMF1 is marginally consistent with observation if 
4 x 10 7 < age/yr < 7 x 10 7 and no more than 30% of 
the ionizing photons are absorbed by dust. Increasing 
/agn much beyond 0.5 is not permitted since there are 
then too many ionizing photons given the bolometric lu- 
minosity and the ionizing photon deficiency problem now 
becomes an ionizing photon excess problem. For example, 
region C on Figure 1 is excluded if IMF1 is valid. 

3.1.2. IMF2 

This steep IMF can only be consistent with observation 
if the age is very young (< 5 x 10 6 yr) and /agn is finely 
tuned to have a value of about 0.5 and if the fraction of 
ionizing photons absorbed by dust is finely tuned to have 
a value of about 0.2 (but not 0). We therefore regard it 
as unlikely. This IMF can never by itself generate a high 
enough value of L^ Q \/M^ to be consistent with observation. 
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Increasing the AGN fraction helps a little but increasing 
it by too much creates an ionizing photon excess problem 
as described in the previous section. 

3.1.3. IMF3 

If we adopt the IMF with a smaller upper-mass limit 
(M u = 3OM : IMF3), the evolutionary loci shift to the 
upper right on Figure 1. In the absence of either absorp- 
tion of ionizing photons by dust or an AGN, the obser- 
vational constraints are satisfied if 6.3 x 10 6 < age/yr < 
1.1 x 10 7 . An AGN fraction greater than zero or an appre- 
ciable fraction of ionizing photons being absorbed by dust 
would generate an ionizing photon excess problem. 

3.2. Current star formation rate in Arp 220 

Since we adopt the constant star formation model, the 
star formation rate (SFR) is coupled uniquely to the bolo- 
metric luminosity for each model. Assuming the SFR of 
1 Mq yr -1 , the ionizing photon production rate is about 
Q = 4.6 x 10 52 photons s" 1 at an age of 1 x 10 7 yr for Model 
with IMF3 (i.e. where it intersects with Box A on Fig. 1). 
This corresponds to a total SFR of the hidden starburst 
of 267M Q yr" 1 . In Arp 220, there are two concentrations 
of thermal dust emission with which the hidden starbursts 
are associated. Downes & Solomon (1998) evaluated the 
luminosity ratio between the western and the eastern nu- 
cleus as 1.5. We assume that the ratio of SFR between 
them is the same as this luminosity ratio. The SFR of 
the western nucleus is then 160M Q yr -1 and that of the 
eastern nucleus is 107M Q yr" 1 . 

The results in this and the previous section depend in- 
tricately on our assumption of a star-formation rate that 
does not vary strongly with time (Point 1 in Section 2). 
An alternative way of explaining the deficiency is to have a 
star-formation rate that was extremely high at some time 
in the past but is low now. The highest mass stars, which 
produce most of the Lyman continuum, would since have 
died (these have very short lifetimes) and therefore the ob- 
served value of Q is low. If the starburst was sufficiently 
powerful to generate many intermediate-mass stars, these 



could produce the observed Lboi- Very large numbers of 
such stars would be required, since each star individually 
contributes only modestly to Lboi- Such a scenario may 
then need to be finely tuned in terms of both the time- 
variation of the star-formation rate and the IMF so as to 
avoid overproducing the total mass in stars (this must be 
less than the dynamical mass), but this is still an impor- 
tant caveat to our results. 

3.3. Visual Extinction 

Once we have chosen a particular model, we can com- 
pute the visual extinction necessary to hide the hidden 
starbursts in the if-band HST image of Scoville et al. 
(1998). We find that Ay must be larger than 30 mag for 
the western nucleus and 40 mag for the eastern nucleus 
for cither IMF1 where more than 0.4 of the ionizing pho- 
tons are absorbed by dust or for IMF3 where none of the 
ionizing photons are absorbed by dust (see the discussions 
in Sections 3.1.1 and 3.1.3). This result is consistent with 
the recent MIR observations; i.e., Ay of 45 mag (Genzel 
et al. 1998). 

4. CONCLUSION 

In all the cases that we consider, the hidden power 
source due to star formation at the 50% level or more in 
order to avoid an ionizing photon excess problem, and this 
starburst must be young (< 7 x 10 7 yr). We derive a cur- 
rent star formation rate of 267M Q yr -1 , and an extinction 
Ay > 30 mag to this hidden starburst. 
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Table 1 
Parameters of IMF 



IMF 


a 


M u 


Mi 


1 


2.35 


100 


1 


2 


3.30 


100 


1 


3 


2.35 


30 


1 




Fig. 1. — The evolutionary locus on the -Lbol/Q — L bo i/M* plane for the IMF models listed in Table 1 from 10 yr (the largest 
L bo i/M* on each locus) to 10° yr (the smallest L bol /M* on each locus). The circles on the line mark intervals of 10 7 yr. The black 
dotted lines show the regions consistent with observation: 8.1 X 10~ 44 LQ(photons s — 1 )~ 1 < L bo i/Q < 1-2 X 10 -43 LQ(photons s -1 )" 1 and 
400-Lq/Mq < L bo i/M* < 600-Lq/Mq. An AGN fraction /agn = is assumed for the black lines. The red and blue lines show the region 
on this diagram consistent with observation for /agn = 0.5 and 0.75 respectively, assuming that the AGN contributes only to L bo i (see the 
text for details). The regions A, B, and C are the regions consistent with observation for /agn values of 0, 0.5, and 0.75. 



